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The properties of pairing correlations in symmetric nuclear matter are studied in the 
rclativistic mean field (RMF) theory with the effective interaction PK1. Considering 
well-known problem that the pairing gap at Fermi surface calculated with RMF effective 
interactions are three times larger than that with Gogny force, an effective factor in the 
particle-particle channel is introduced. For the RMF calculation with PK1, an effective 
factor 0.76 give a maximum pairing gap 3.2 MeV at Fermi momentum 0.9 fm , which 
are consistent with the result with Gogny force. 



1. Introduction 

The mean field theory, including non-relativistic mean field theory with effective 
nucleon-nucleon interactions such as Skyrme or Gogny, and the relativistic mean 
field (RMF) theory, has received wide attention due to its successful descriptions 
of lots of nuclear phenomena during the past years. In the framework of the RMF 
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theoryP the nucleons interact via the exchanges of mesons and photons. The rep- 
resentations with large scalar and vector fields in nuclei provide simpler and more 
efficient descriptions than non-relativistic approaches that hide these scales. In these 
sense, the RMF theory is more fundamental. With a limited number of free parame- 
ters, i.e. meson masses and meson-nucleon coupling constants, the RMF theory has 
proved to be successful in quantitatively describing the prop erties of nuclear matter 
and neutron starsP nuclei near the valley of stabilitypHEl anc j exotic nuclei with 
large neutron or prot on exce ss with proper treatment of the pairing correlations 
and continuum effects P^^^ It giv es naturally the spin-orbit potential, th e origin 
of the pseudo-spin symmetry^SLI] as a relativistic symmetrypHHHIHl an( j t ne 
spin symmetry in the anti-nucleon spectrumP^ 

Since 1950's, a large number of striking experimental facts, such as the binding 
energy difference between even-even and odd-even nuclei and a systemic reduction 
of the moments of inertia of even-even nuclei compared with their neighboring odd- 
even nuclei in def ormed nuclei, suggest the existence of the superfhud phenomena in 

FTT1 • • m^i 

these systems. 1 ^ In astrophysics, the origin of pulsar glitches, 1 ^ which are sudden 

discontinuities in the spin-down of pulsars, can also be understood via the super- 
fluidity in the inner crust of these stars. All of these phenomena suggest that the 
pairing correlations play an important role in theoretical study on the properties 
of nuclei structure and nuclear matter. 

The first relativistic study of superfluidity in infinite nuclear matter was done by 
Kucharek and Ring in 199ll^They studied the pairing correlations in symmetric 
nuclear matter using the relativistic Hartree-Bogoliubov (RHB) method with an 
one-boson-exchange (OBE) potential in the particle-particle channel. However the 
resulting pairing gap at t he Fermi surface are about three times larger than that 
with the Gogny forcePS Therefore, the effective pairing interaction used in RHB 
calculation are either the finite range Gogny force or the Skyrme type zero-range 
force. In particular, with the Skyrme type zero-range force, lots of achievements 
have been made to describe the nuclei far away from the line of /^-stabi lity with 
proper treatment of the pairing correlations and continuum effects P l ^^ 

However, it is still an open problem how the same nucleons interaction via 
the exchanges of mesons and photons in the Hartree channel can be used in the 
particle-particle channel as well. In fact, the large pairing gap in RHB calculation 
with an OBE potential in the particle-particle channel comes from the behavior 
of the pairing matrix elements at large momentaPH The various effective forces 
in RMF models are adjusted for mean-field calculations in Hartree channel only, 
i.e., only for momenta below the Fermi momentum, thus a realistic particle-particle 
interaction can have very different behavior at high momenta. Therefore, in order 
to get reasonable values for the pairing gap, one can use a suitable value of cut off 
in the momentum space in the relativistic mean field calculationspEIESES] or 
consider various effects, such as the medium polarization, the in-medium meson 
mass decrease, and the mesons nonlinear terms to reduce the pairing gap in nuclear 
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In this paper, the properties of pairing correlations in symmetric nuclear matter 
will be studied in the RMF theory with the newly developed effective interaction 
PKlP3 In order to solve the well-known problem that the pairing gap at Fermi 
surface calculated with RMF effective interactions are three times larger than that 
with Gogny force, an effective factor in the particle-particle channel will be intro- 
duced. In section II, a brief description of the RMF theory and RHB theory for 
nuclear matter is presented. The results and discussions are given in the following 
section. Finally in the last section, a brief summary is given. 



2. Theoretical Framework 

2.1. Relativistic Mean Field Theory 

A general review of the RMF theory and its application in nuclear physics could 
be found in Refs. [21 [H [5]. Here a brief review of RMF theory for nuclear matter is 
given. The start point of RMF theory is an effective Lagrangian density with the 
nucleons interacting via the exchange of various mesons and the photon: 



il^d^ - to - g a a - g^u^ - g p ^T ■ p p - 



2 p 



1 



+ ^ad^a - ^m 2 a a 2 - U{a) - ^ v % v + ^w"^ + U(uj) 



(1) 



The Dirac spinor ip denotes the nucleon with mass to. The isoscalar scalar a- and 
isoscalar vector u-mesons offer medium-range attractive and short-range repulsive 
interactions respectively, and the isovector vector p-meson provides the necessary 
isospin asymmetry. Their masses are denoted by m a , to w and m p . g a , g u and g p are 
the corresponding meson-nucleon coupling constants, r is the isospin of nucleon, 
and T3 is its three-component. The nonlinear a and ui self-interactions U (a) and 
U(u>) are respectively denoted as, 



1 



1 



U(o) = -g2<J A + -g 3 cr\ U(w) = -c 3 K^, 



(2) 



with the self-coupling constants <?2, ff3 and C3. The field tensors fl^, R pv , and A pu 
are as the followings, 



d p uo„ - d u uo p , R pu = dfip u - dvp^, A^ = d^A^ - d v A h 



(3) 



The classical variation principle gives the following equations of motion, 

1 - To 1 

i^dn - to - g a a - g^u; p - g p ^f ■ p p - - A p tp = 0, (4) 
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for the nucleoli spinors and 

{d»d p + ml)a = -gM ~ g^ 2 - g^\ (5) 
d^ v + mlu v = g u ftfl> ~ cstrfuf) 3 , (6) 
d p R^ + my = g p tfiS> + g P P P x ST, (7) 

8 P A^ = efrfl^p-^ (8) 

for the mesons and photon, where the sum over all the particle states in the no-sea 
approximation is adopted for the source term in the Eqs. ([51 El El EJ)- 



2.2. Relativistic Hartree-Bogoliubov Theory 

Usually, in the RMF theory, the mesons are treated as classical fields. In order to 
describe the superfluidity of the nuclear many-body system, one needs to quantize 
not only the nucleon but also the meson fields. By using the well-known canonical 
quantization method and the Green's function techniques, neglecting retardation 
effects and the Fock term as it is mostly done in RMF, one can get the so-called 
Relativistic Hartree-Bogoliubov (RHB) equation}!^ 

where 

h= a-p + V + /3(M + S), (10) 
is the Dirac Hamiltonian with the scalar potential S and vector potential V, 

S = g a a, V = p + g p ^r ■ p p + ef^p-A^ . (11) 

The pairing field is 

A-ab — ~ ^2 VabcdHcd, (12) 
cd 

where V a bcd is the two-body effective interaction in the particle-particle (pp) chan- 
nel and the pairing tensor n a b = Y^k^ak^bk- The quasi-particle eigenvectors are 
denoted as 14), and is its corresponding quasi-particle energies. 

The chemical potential A in Eq. §§§ is determined by the particle number with 
the subsidiary condition, ^ fe V fc 2 = N. 



2.3. Application to Symmetric Nuclear Matter 

For the static, uniform infinite nuclear matter, the coulomb field is neglected, and 
the space-like components as well as the differential of the time-like components of 
the meson fields vanish. Furthermore, for symmetric nuclear matter, the p meson 
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has no contribution on the mean field potential. Then the scalar potential S and 
vector potential V are constants and have the simple form, 

S = g a (a), V = g u (u> ). (13) 

The RHB equation © can be decomposed into (2 x 2) matrices of BCS-typeP^ 

e(k)-X A(fc) \(u(k)\_ e(k) (u(k)\ 
A(fc) -e(fc) + A J \v(k) J ~ () \v(k) ) ' 1 ' 

where the eigenvalue of the Dirac Hamiltonian for positive energies is denoted as, 
e(k) = V + ^k 2 + (m + S) 2 , the Fermi energy A = e(kp), and the quasi-particle 
energy e(fc) = U (e(k) — A) 2 + A 2 (fc). The corresponding occupation numbers v 2 (k) 
have the form 

Ak) = l M- , £{k) ~ x ). (15) 



The pairing field A(fc) obey the usual gap equation, 

where / is an effective factor introduced to reduce the pairing potential. The effec- 
tive interaction in the pp channel v pp (k,p) is the one-meson exchange potential, 

v PP {k,p) = Vp p (k,p) + v pp {k,p) + vP p (k,p), (17) 

where 

(j / ,s gg f ( £ *(p)^ £ *(fc)) 2 +m 2 -4m* 2 (fc+p) 2 + m 2 1 



3 5 2 £ *(fc) £ *(p)-m* 2 1 (fc+p) 2 



m: 



v u (v k) = — v ' — — - ln v ' u (19) 

"ppKP'K) £ *( fc ) £ * (:p ) 2 pk (k-p) 2 + ml' [ ' 



vP(pk) - 9 * 2e*(k)e*(p)-m* 2 (k+ P ) ...„ 
u pp [p ^> £ *(k)e*{p) 2pk {k-p) 2 + mf K 1 

with effective mass to* = m + .9<xC, and e*(fc) = \/k 2 + to* 2 . 

The meson fields are replaced by their mean values, and could be solved from 
the corresponding equations of motion by the various giving nucleon densities, 

m l<* = -g*Ps - 92V 2 - (21) 
to 2 w = 9uPv - c 3 Wo, (22) 

where p s and p v are respectively the scalar- and baryon-density, 

2 r°° m + g a a ..2 fu \u2. 
lo v/fc 2 + (m + g a a) 2 

Pv = ip^ip = — / v 2 (k)k 2 dk. (24) 



Ps = H> = — / , ?2 == =v\k)k'dK (23) 
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3. Results and discussions 

For a given Fermi momentum Uf, the coupled Eqs. (I14 [[T61[2T|[22|) can be solved self- 
consistently by iteration. The properties of 1 S$ pairing correlations of symmetric 
nuclear matter are studied with the newly developed effective interaction PK1 which 
takes into account the self-interactions of cr-meson and w-meson as well as the 
isospin dependence of the nuclear matter.^ 

The momentum integration in the gap equation, in principle, should go to in- 
finity. In actual calculations, it is necessary to have a cutoff in the momentum 
space and the convergence of the pairing gap on the cutoff momentum should be 
checked. The dependence of the pairing gap on the cutoff momentum for different 
Fermi momentum is given in Fig. [1] The results indicate that the cutoff momentum 
Kc > 10 fm _1 will guarantee the numerical convergence. In the following, Kc = 20 
fm will be adopted and the corresponding effective interaction in the pp chan- 
nel, the momentum-dependence of the pairing gap, and the influence of effective 
interactions on the pairing gap at the Fermi surface, etc., will be investigated. 




Fig. 1. The pairing gap A(kp) at the Fermi surface as a function of the cutoff momentum kc in 
the momentum space for different values of Fermi momentum kp = 0.3, 0.6, 0.9, and 1.2 fm 
with the effective interaction PK1. 



3.1. The effective interaction in the pp channel 

The contour plot for the effective interaction in the pp channel v pp (k,p) for different 
Fermi momentum with the effective interaction PK1 are shown in Fig. [21 where the 
contours with negative values are denoted by red dashed lines. 
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The interaction is attractive for small and repulsive for larger momenta k and 
p, or equivalently, attractive for large distances and repulsive for small distances. 
At around 1.5 fm -1 , the interaction will change from being attractive to being 
repulsive. The repulsive interaction reaches its maximum value at the momenta k 
and p around 4 fm -1 . The maximum repulsive interaction increases with the Fermi 
momentum. 

The behavior of the effective interaction in the pp channel v pp (k,p) can be un- 
derstood from contributions of different mesons, as shown in Fig. [3] The scalar 
meson a provides the attractive part of the effective interaction, v pp (p,k), with a 
peak value at zero momentum, and approaching zero with the increasing momen- 
tum. While the vector meson lo and p provide the repulsive part, which extend to 
higher momenta than v pp {p, k). The main contribution for the repulsive part comes 
from the to meson, as seen in the figures, v pp (p, k) is one order of magnitude larger 
than vP p (p, k). 

Different mesons contributions to the effective interaction in the pp channel 
v pp (k,p) as a function of p at k = &f = 0.9 fm -1 with the effective interaction PK1 
are shown in Fig. |U The sum of all the mesons contributions has considerable re- 
pulsive contributions for momenta larger than about 1.5 fm -1 , which is remarkably 
different from the one by the Gogny force calculation}!^ 
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Fig. 2. Contour plots for the effective interaction in the pp channel v pp (k,p) as a function of the 
momenta p and k for different Fermi momentum kp = 0.3, 0.6, 0.9, and 1.2 fm — 1 with the effective 
interaction PK1. The contour lines have a distance of 1 (fm 2 ) and the negative values are denoted 
as red dashed lines. 
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Fig. 3. Contour lines of different mesons contributions for the effective interaction in the pp channel 
v pp (k,p) as a function of the momenta p for Fermi momentum kp = 0.9 fm — 1 with the effective 
interaction PK1. The contour lines have unit of fm 2 and the negative contours are denoted as red 
colors. 




Fig. 4. Different mesons contributions to the effective interaction in the pp channel v pp (k,p) as 
a function of p at k = kp = 0.9 fm — 1 with the effective interaction PK1. The dashed lines 
corresponds to that of the a-, u>- and p-fields, the solid line represents the total contribution. 
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3.2. The momentum- dependence of the pairing gap 

For a given v pp (k,p), the momentum-dependence of the pairing gap A(/c) can be 
obtained from Eq. (|16p. The momentum-dependence of A(fe) with different Fermi 
momentum is shown in Fig. [5l The pairing gap has large and positive values at 
small momenta, then decreases with the momentum and changes its sign around 
2.5 fm -1 . It continues to decrease till ~ 4.0 fm -1 , then slowly turns to zero. The 
difference of pairing gap for different Fermi momentum is revealed mainly by the 
behaviors at low momentum. At zero momentum, the largest pairing gap occurs at 
the Fermi momentum 0.9 fm -1 . 




Fig. 5. The pairing gap A(fc) as a function of the momentum k for Fermi momentum fcj?=0.3, 0.6, 
0.9, and 1.2 fm - 1 with the effective interaction PK1. 



3.3. The pairing gap at the Fermi surface 

One of most important properties of pairing gap is its value at the Fermi surface. 
In Fig. [HI the pairing gap A(kp) at the Fermi surface as a function of the Fermi 
momentum kp with the effective interaction PK1 is shown in comparison with the 
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results obtained with the effective interactions NLlpHl NL2,^21 NL3J321 NLSHpS 
TMllS^! and the results calculated with Gogny force and Bonn potential!^ 

It is found that the pairing gap A(fc^) is strongly dependent on the nuclear 
matter density, or equivalently, the Fermi momentum. The pairing gap A(/c_f) in- 
creases with the Fermi momentum (or density) and reaches maximum at Fermi 
momentum fci? ss 0.9 fm , then rapidly drops down to zero. Usually, the pairing 
gap at the Fermi surface calculated with RMF effective interactions is almost three 
times larger than the value calculated with Gogny force. Moreover, the pairing gap 
at lower Fermi momentum does not vanish in calculations with the RMF effec- 
tive interactions, while it docs vanish in calculations with Gogny force or Bonn 
potential. 

These differences come from the integral in Eq. (|16|) for the pairing gap, which 
depends on the products of pairing gap parameter A(p) and the effective interaction 
in the pp channel v pp (k,p). From v pp {k,p) in Fig. [2] and A(p) in Fig. [5j it is found 
that considerable contributions to the integral in Eq. (|16p may come from high 
momenta region. While the various effective forces in RMF models are adjusted 
for mean-field calculations in Hartre e ch annel only, i.e., they are only valid for 
momenta below the Fermi momentunji^ and the a realistic interaction in the pp 
channel v pp (k,p) can have very different behavior at high momenta. 

Considering that the RMF effective interactions give a too much strong pairing 
field, an effective factor is introduced in the particle-particle channel to reduce the 
pairing gap. For PK1, if a factor / = 0.76 is introduced, the resulting pairing gap 
is almost the same as those with Gogny force or Bonn potential, and a maximum 
pairing gap 3.2 MeV is obtained at Fermi momentum 0.9 fin , as shown in Fig. [6j 
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Fig. 6. The pairing gap A(kp) at the Fermi surface as a function of the Fermi momentum kp for 
different effective interactions. The results of Gogny and Bonn come from the Ref. [21]. 
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4. Summary 

The pairing properties in the 1 S'o channel for symmetric nuclear matter have been 
studied in the RMF theory with the effective interaction PK1. The one-meson 
exchange potential is used in the particle-particle channel in consistent with the 
particle-hole channel. The effective interaction in the pp channel is found to be 
attractive at small momenta with attractive range around 1.5 fin -1 and repulsive 
at large momenta. The pairing gap at the Fermi surface is strongly dependent on 
the nuclear matter density. It grows as the Fermi momentum increases, reaching its 
maximum values at Fermi momentum around 0.9 fm -1 , then drops down to zero 
rapidly. Considering the fact that the pairing gap at Fermi momentum calculated 
with RMF effective interactions are three times larger than that with Gogny force, 
an effective factor in the particle-particle channel is introduced. For the effective 
interaction PK1, a factor / = 0.76 will produce almost the same results as those 
with Gogny force or with Bonn potential. 
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